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^ : Abstract 

This paper investigates the possible use of the Hyperspherical Adiabatic basis in the description 
of scattering states of a three-body system. In particular, we analyze a 1+2 collision process below 
the three-body breakup. The convergence patterns for the observables of interest are analyzed by 
comparison to a unitary equivalent Hyperspherical Harmonic expansion. Furthermore, we compare 
and discuss two different possible choices for describing the asymptotic configurations of the system, 
^^ . related to the use of Jacobi or hyperspherical coordinates. In order to illustrate the difficulties and 



advantages of the approach two simple numerical applications are shown in the case of neutron- 



C^ ' deuteron scattering at low energies using s-wave interactions. We found that the optimization 

0|0 , driven by the Hyperspherical Adiabatic basis is not as efficient for scattering states as in bound 



state applications. 



I. INTRODUCTION 

The Hyperspherical Adiabatic (HA) method is based on the parametrization of the inter- 
nal degrees of freedom with hyperspherical coordinates (see Refs.|l] and references therein). 
The method then consists in expanding the system's wavefunction on a basis made of hyper- 
angular optimized functions (the adiabatic basis set) times (unknown) hyperradial functions. 
The hyperangular basis elements are taken as the eigenvectors of the Hamiltonian operator 
for a fixed value of the hyperradius p. Once those eigenvectors have been calculated, the 
hyperradial functions are obtained as the solutions of a system of coupled one-dimensional 
differential equations. The advantages of such approach are that the HA basis should drive 
a quick convergence for the expansion, due to its optimization, the payback is represented 
by the necessity and the difficulty in calculating accurately the first and second derivatives 
of the adiabatic basis set with respect to the hyperradius. Those terms are crucial to the 
method as they represent the coupling terms between the various hyperradial differential 
equations. In some applications of the HA method it was shown that the strong coupling 
between pair of elements of the adiabatic basis makes the hyperradial problem particularly 
hard to solve [2(]. 

The properties of the adiabatic basis functions have been object of several studies and 
are well-known. In particular, in the asymptotic limit of large hyperradius the HA functions 
are known to converge towards the scattering states of the three-body system, both below 
and above break-up. This characteristic makes the adiabatic expansion a valid choice to 
describe the three-body continuum states. In the literature there are several studies of the 
bound spectrum of a three-nucleon system by means of the HA method |3|, y, |5| , but very 
few dealing with continuum states [6]. This paper investigates the possibility of using the 
HA approach to describe a three-body elastic process in which a particle collides the other 
two, initially forming a bound state. The object of this work is the study of the appropriate 
boundary conditions to be imposed to the hyperradial functions as p ^ oo and a careful 
analysis of the convergence properties of the HA expansion. 

In order to quantitatively understand the pattern of convergence of the HA expansion we 
make use of the parallelism that can be built between the HA method and the Hyperspherical 
Harmonic (HH) expansion. In fact, we can consider two different expansions for the system's 
wavefunction, one in terms of Na HA basis elements, and the second in terms of Nh HH basis 



elements. When Na = Nh the two expansions are connected by a unitary transformation 
and therefore must yield identical results. Since the HH basis has been used several times to 



3,3 



describe scattering states [7|, |8| , we exploit this knowledge to study the convergence of the 
HA expansion. In particular, we will study the convergence properties of the L = phase 
shift at low energies in a 1 + 2 collision, which has been used as a benchmark problem in 
literature (see for example J9|). 

The problem of the boundary conditions to be imposed to the hyperradial functions is 
related to the difficulties associated with obtaining the eigenvectors and eigenvalues of the 
adiabatic Hamiltonian at large values of the hyperradius. As the lowest adiabatic functions 
tend to the two-body bound wavefunctions, an accurate description of those states using, 
for example, the expansion in HH functions is known to be very difficult. This is because, 
as p — > oo, the two-body bound states are localized in a very small zone of the hyperangular 
phase-space. Consequently, this particular configuration necessitates a large number of HH 
functions to be described 21]. In fact, it can be shown that the number of HH required to 
reproduce this type of spatial configuration grows exponentially with the hyperradius. If 
the interest is limited to study deep three-body bound states, the problem just described 
does not manifest and a tractable number of HH functions suffices for a good accuracy. 
Due to the finite hyperradial size of the associated wavefunction, the adiabatic Hamiltonian 
needs to be solved only up to a non so large value of the hyperradius. However, there 
are cases in which shallow bound states are present (as Efimov states) and the adiabatic 
Hamiltonian needs to be solved for very large values of p. Furthermore, for energies in the 
continuum, the associated three-body scattering wavefunction has an infinite extension and 
a direct application of the HA necessitates of the solution of the adiabatic Hamiltonian at 
very large values of p, too. In order to obtain accurate asymptotic solutions to the adiabatic 
Hamiltonian we have followed in detail the procedure outlined by Nielsen and co-workers [l| . 

Finally, interest in this work is also sparkled by an article of Fabre de la Ripelle Gj , where 
he suggested the possibility of expanding the three-body asymptotic scattering states into 
the adiabatic basis set, and retaining only the first term in such an expansion, resulting in 
a considerable reduction of the numerical burden. We will analyze this truncation together 
with the contribution of higher terms. 

This paper is organized as follows: in the next section the HA method is presented, by 
first introducing the notation. The expansion of the HA basis in terms of the HH functions 



is given as well as the method to describe the HA functions and the adiabatic potentials 
at large values of p. Section III treats the problem of scattering states. Two different 
methods of implementing the Kohn Variational Principle are given in conjunction with the 
HA basis. The asymptotic conditions are given in terms of the distance between the incident 
particle and the two-body system and in terms of p. Section IV is devoted to numerical 
applications. Results are presented using a simple Gaussian two-body potential and the 
semi-phenomenological s-wave MT-III potential [lOj. The final Section is devoted to the 
conclusions and perspectives. 

II. HYPERSPHERICAL ADIABATIC METHOD 

Let us consider a system of three identical particles of mass m, in a state of total orbital 
angular momentum L = 0. Other quantum numbers are represented by the total spin S", 
total isospin T, and the symmetry under particle permutation H, which can take the values a 
(anti- symmetric, for three fermions) or s (symmetric, in the case of three bosons). A further 
quantum number needed to uniquely identify each wavefunction is given by the vibrational 
number n (n = 1, 2, ...) for bound states or the energy E for continuum states. 

Let us start from the definition of Jacobi coordinates {xj,yj} 



Xi = —^[J^j - TfcJ 



1 

Yi = -T^i^j + ^k - 2ri) (1) 

where {rj} are the Cartesian coordinates of the three particles and i,j,k = 1,2,3 cyclic. 
The hyperspherical variables {p, 6i} are defined as follows 

Xi = p cos 9i, Yi = p sin 9i (2) 

where p is the hyperradius which is symmetric under any permutation of the three particles 
and 6i is the hyperangle, which is dependent on the particular choice of the Jacobi coordinate 
system. In terms of the interparticle distances r^ = |rj — rj| = y/2Xk the hyperradius reads: 



P=-7^\/rh + r^3 + rli ■ (3) 

In addition to p and 6i there are four more coordinates needed to parametrize all the possible 
spacial configurations of the three particles, for example the four polar angles which define 



the orientation of the two Jacobi vectors with respect to the laboratory frame of reference. 
However, in the particular case of total orbital angular momentum L = 0, the number of 
such coordinates can be reduced to just one non-trivial functional dependence, represented 
by the cosine yUj of the angle between the two Jacobi vectors {xj,yj}: 

/ii = Xi ■ yi/{XiYi}. (4) 

In the following we will refer to the set of hyperangles {6'j,/ij} as Qi, or more in general as 
Q = {6, yu} when there is no need to specify the choice of a particular permutation of the 
particles defining a set of Jacobi coordinates. 

The Hamiltonian operator Ti takes the following expression in hyperspherical coordinates 

where V is the potential energy operator, Tp is the hyperradial operator 

^ <P b d . , 

and G^ is the grand-angular operator 

Ad d f^ P 

VT^^dz^ ' dz^ cos^6^ sm'^Q- ^^ 

where z = cos 26 and £x and iy are the angular momentum operators associated with the x 



and y vectors, respectively. The volume element is p'dp\J\ — z'^dzdjj,. 

The system wavefunction \1/, with quantum numbers L, S, T, U, and n (or E), is expanded 
as follows: 

oo 

where {$^'^^^} are the eigenfunctions of the operator Tin made of the hyperangular part of 
the kinetic operator plus the potential energy operator, in which p acts only as a parameter: 

Hn$^^^" = [^G^ + V^j $^"^"(P, n) = t/.(p)<|.^^n(p, n). (9) 

The set of eigenfunctions {$^'^'^} is known as the adiabatic basis set, and the associated 
eigenvalues {U^j^p)} as the adiabatic curves or potentials. In practical calculations, the 
infinite expansion of eq. (IHl) needs to be truncated to a finite number of basis elements, say 
A^^. The convergence for the observables of interest with respect to this parameter is then 
checked. 



The initial Hamiltonian problem is thus tackled in two steps: firstly, the HA basis func- 
tions {^^^^} and the associated potentials {Uu{p)} are calculated by solving eq. Q. Sec- 
ondly, the hyperradial functions m"(p) are obtained as the solutions of a system of A^^ 
coupled one-dimensional differential equations, which can be expressed as follows ll[ |: 



Na 

E 

u=l 






dp 



d 



dp 



{u' = 1,...,Na), (10) 

where the coupling terms B,yi^, C^ii, follow from the dependence on p of the HA basis : 

n^ d^,, d^ 

and 



Cu'u{p) 






mp'^ 



dp 



(12) 



For bound states solutions, and short range potentials, the functions {uy} tend to zero 
exponentially as p ^ cxd, whereas for scattering states the boundary conditions to be imposed 
to the {uy} will be discussed in the next Section. 

The first step in the implementation of an HA calculation consists in obtaining the adi- 
abatic basis elements and the associated adiabatic potentials, solutions of eq. ([9]), for a 
number of values of p. Among several available techniques we have chosen to use a varia- 
tional approach, by expanding the functions {<l'^'^^^} onto a set of Hyperspherical Harmonics 
(HH) of size Nh. In order to define a basis set with the desired properties under particle 
permutation, we combine opportunely hyperspherical polynomials based on different set of 
Jacobi coordinates J8|. The expansion for $^•5'^^ j-g^^jg; 

Nh 

J2DUp)\kl,LSTIl), (13) 

kl 
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LSTU 



with the basis element given, for L = 0, by 

3 

\kl, OSU) = J2 l^^'^H^^i) ^Ti^Si 



(14) 



1=1 



where Si (Tj) indicates the coupling of particles jki to a state of total spin 5* (total isospin 
T), and the hyperspherical polynomial is written as (see for instance Ref. 12, ll3|] for more 
details): 

(^)p^''(fi) = N,,{1 - z'f/'^Pi^'/'^'^'/\z)Pi{p) , (15) 



where P^ is a Jacobi polynomial, Pi is a Legendre polynomial and Nki is a normalization 
factor. The HH so defined are eigenfunctions of the grand-angular operator, 

G^|A;/,O^Tn) = K{K + A)\kl,OSTU), (16) 

where K is the grand-angular quantum number {K = 2k + 21) . 

The unknown coefficients {-D^/} in eq. (1131) . and the adiabatic potential {U^} are obtained 
as the eigenvectors and eigenvalues, respectively, of the following generalized eigenvalue 
problem 

Nh 

^{k'l', LSTUlHn - U\kl, LSTU)Dki = 0. (17) 

kl 

In practical calculations the size A'^^ of the HH basis set is increased until convergence 
is reached for the desired number A^^ of adiabatic potentials {U,y}. However, it is well 
known that the convergence becomes harder to achieve the larger the value of p. The reason 
for this behavior is connected to the specific properties of the HA basis set at large p. 
Namely, the lowest adiabatic potentials tend to the binding energies of all possible two- 
body subsystems, and the associated HA basis elements to the two-body wavef unctions, 
opportunely normalized. The HH expansion is not optimal for reproducing wavefunctions 
with similar characteristics, which become the more localized the larger p . This convergence 
problem can be further enhanced by the presence of a hard core repulsion in the two-body 
potential. If the calculation is to be limited to the three-body bound states, and in absence 
of very extended ones such as the Efimov states, the limited radius of convergence of the 
HH expansion does not constitute a problem. When the calculation is extended to the 
continuum energy region, however, the accurate determination of the adiabatic curves and 
functions at very large p becomes essential for the convergence of the results. In order to 
overcome this problem Blume and co-workers ^] advocate the use of splines, which at large 
p converge significantly faster than the HH. Alternatively, when p is much larger than the 
range of the two-body interaction, approximations for the HA basis elements and potentials 
can be obtained by solving a non-homogeneous one- dimensional differential equation. A 
brief illustration of this second approach is summarized below, based on the work of Nielsen 







and co-workers [l|. Let us start from the definition of the reduced amplitudes 
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Z^ u Z^cos^,sin^/ ^ ^ 



=1,3 i=l 



each one having the set of quantum numbers LSTU. They are the solutions of the Faddeev 
equations, that for s-wave potentials read 



+ Vi^pcosOi) - K{p) (l),{p,9i 



2mp2 d9f 
-cos Oi sin 9iViV2p cos 9i) f dp, ($[/) + <^i^^) (19) 

where Xu{p) = Uu{p) — Ah^ /{2mp'^). Defining tq = -\/2pcos^o the range of the (short- 
range) potential, we observe that, for large values of p, the potential V{^/2p cos 9i) can be 
considered different for zero only for values of 9i in the interval ^o ^ 0i ^ 7r/2, which is 
the smaller the larger p. Accordingly, the above equation has two regimes depending the 
values of 9i. It is homogeneous for 9, < 9q. For values in which the potential is not zero we 
have to evaluate the non- homogeneous term which depends on the amplitudes j, k. From 
the relation between the different sets of Jacobi coordinates, the region of values of 9i where 
V is different from zero correspond to the values 9j ~ tt/6 and 9k ~ 7t/6. In this region each 
of these amplitudes is governed by the corresponding homogeneous Faddeev equation. For 
example, for the j-amplitude, the possible solutions depending on the value of Ay are 

(j)^{p, 9j) = A sm{K9j) X^>0 
Mp, 0,) = A{e^^'^ - e-^-'^O A, < , (20) 

and similarly for the /c-amplitude, where k^ = 2m\X^\/h?. Replacing these expressions in 
the Faddeev equation (IT^ . its asymptotic form can be obtained: 

(-^^ + nV2pcos0) - A.(p)') Mp,0) = V{V2pcos9)Af{p,9) (21) 

When the equation describes a two-body bound state with a third particle far away, 
Xi, is negative and tends to the two-body bound state energy. The corresponding non- 
homogeneous term is 

k{TT/2-e) _ -fc(7r/2-6») fe7r/6 _ „-fc7r/6 

For positive values of Ajy the adiabatic functions describe asymptotically three free particles 
and 

fip,9) = -'-^^^^sin[k{n/2-9)]/k. (23) 



y4 is a normalization constant to be determined from the solutions. The boundary conditions 
for the functions 0y are 0t.(p, 0) = 0y(p, 7r/2) = 0, which determine completely the solutions 
ofeq. ([2ID. 

In practical applications the adiabatic potentials {U,^} and the HA basis elements {^^^^^} 
are obtained as solutions of eq. (TT71) for p < po and of eq. I^I\i for p > po, respectively. The 
matching point po needs to be chosen larger than the range ro of the two-body potential V. 
There is a zone around the matching point in which, for a sufficient large value of Nh, the 
solutions obtained from the HH expansion or by solving eq. (!2T!) for each value of u become 
indistinguishable from each other. In this way we link the definitions of po and Nh as the 
values for which the solutions of eqs. fll7p and i^T^ can be accurately matched. In fact, if 
the functions (pu obtained by solving eq. fl^Tl) are themselves expanded into the HH basis, 
the coefficients of this expansion can be individually matched to the equivalent coefficients 
obtained through solving eq. (fT7|) for the same value of p. 

In the following we discuss the solutions of the the system of coupled differential equations 
( TTOl) in the case of bound states. The hyperradial functions {mJJ} can be expanded into 



normalized generalized Laguerre polynomials times and exponential function 14l |: 



"(P) = E ^w^^n/?P) exp [-Pp/2], (24) 



m=0 



where /3 is a non-linear parameter which can be used to improve the convergence of the ex- 
pansion [iSij. The coefficients {A'^^^} can be found by means of the Rayleigh-Ritz variational 
principle, whose implementation requires the solution of the following eigenvalue problem: 

"^{m'u'in - E\mu)Amu = 0, (25) 

mu 

where the ortonormalized basis element Imu) is defined as 

\mu) = Lt\Pp) exp [-/5p/2]<|.f ^"(p, fi). (26) 

The size of the variational problem is M = A'p x A^^, where Na is the number of adiabatic 
basis functions retained in expansion of eq. ([H]), and Np is the number of Laguerre polynomi- 
als used in expansion of eq. fl24p . For sake of simplicity all functions m" are expanded using 
the same number of Laguerre polynomials, although this is not strictly necessary. The eigen- 
values {En } {n = 1, 2, . . .) represent upper bounds to the eigenvalues of the Hamiltonian 



problem ([5]) and converge towards them monotonically as M is increased. The associated 
set of coefficients {A^^} provide approximations to the system wavef unctions. 

As it has been mentioned before, there is a complete equivalence between the two methods 
if they include the same number of HH functions. In fact the expansion for ^ in eq. ([H]) can 
be written also as: 

N„ 

vi/f^" = 5^<,(p)|fc/,L^Tn), (27) 



kl 

and from eq. ( TT3l) the following relation can be obtained 

Na 

wUp) = Y,ul{p)Dl,{p) . (28) 

If A^^ is set equal to N^ the matrix Z)^; represents a unitary transformation between the HA 
and HH basis sets, therefore the two expansions must produce identical sets of eigenvalues 
and eigenvectors. Consequently, if in a specific problem, the desired accuracy is reached using 
Nh HH basis functions, the use of a larger number of HH basis elements in the expansion 
of the adiabatic basis functions is superfluous. However, we can expect that the number of 
adiabatic functions Na needed to reach the same accuracy will be A^a ^ Nh- This is because 
the HA functions have been optimized to the specific Hamiltonian problem by solving eq. 
([9]) for each value of the hyperradius. We would like to stress the fact that the equivalence 
between the HH and the HA method using a tractable number Nh of HH functions applies 
in presence of deep bound states. When shallow bound states, as Efimov states, are present 
the situation changes and a direct application of the HH method encounter the problem of 
the inclusion of a very large number of basis states in the expansion of the wavefunction. 
This is related to the correct description of the adiabatic potentials in the asymptotic regime. 
In this case the use of the asymptotic form of the Faddeev equations given above proves to 
be extremely useful, as for example in the solution of three Helium atoms system [l|. 

III. SCATTERING OBSERVABLE CALCULATIONS 

In this section we apply the HA expansion to the study of continuum states of a three- 
body system. The case considered will be the scattering of one particle colliding other two 
forming a dimer, at energies below the three-body breakup threshold. The wavefunction for 
the system can be written as 

^ = ^c + ^a, (29) 

10 



where the first term is C? and describes the system configurations in which the three particles 
are all close to each other. The second term represents the solution of the Schroedinger 
equation in the asymptotic region in which the incident particle does not interact with the 
other two ( the discussion will be limited to short range potentials). Moreover, we will 
consider the case of a two-body interaction that supports only one dimer bound state of 
energy E"^^ . Accordingly, we will consider energies E"^^ < E < 0. 

The explicit form of the term \l'a depends on the energy E of the system. However, the 
particular choice of the function \l'a is rather arbitrary, as it can be modified by adding or 
subtracting any C^ function. In the following we will consider and compare two different 
expressions for the asymptotic function \t'a. Practical applications will be shown for the 
case of nucleon-deuteron scattering using the semi-realistic s-wave MT-111 potential, as the 
repulsive core of the potential allows a better understanding of the numerical problems 
associated with the method's implementation. 

A. Scattering below Break-up: Method 1 

The \E'a term must describe the asymptotic state of the dimer plus a third particle. 
Therefore, the most natural choice for this term leads to building two independent and 
symmetrized states, that for L = 0, read as follows: 



.i? _ '^ X r 9 (rj) sin [kyV i 



^ST = E^^^^nr^^o(/i.)|5T), (30) 



and 

o^ ^^ xfA ( ^ cos[/i:^t/i](l-exp[-7t/^]) 

. l^yVi 

The distance between particle i and particles j, k forming a dimer is |/j, 0rf(r) is the dimer 
wavefunction of energy E"^^ , k'l = 4m{E — \E'^^\)/3h'^ and A/" is a normalization factor chosen 
so that 

\ilorp\rL E\ilQrp/ (ilorplri E \il grpj = 1/2 . ("2) 

The behavior of the function Qgj^ for yi -^ has been regularized by means of an opportune 
factor. The constant 7 can be consider a non linear parameter of the scattering wave 
function. The final result should be independent of the value chosen for it but a wrong choice 
can slow down the convergence significantly. A reasonable choice could be 7 ?^ y^m\E'^^\/h?. 

11 



A general scattering state is given by defining the following linear combinations 

^ST = UqR^sT + Uoi^sT^ (33) 

and 

^ST = uirQst + uuQsj^. (34) 

The term "^a, having total spin S and total isospin T, can thus be written as 

■^a = ^°ST + ^^It (35) 



where different choices for the matrix u can be used to define the scattering matrix C 16|. 
Here we will use 

defining C = S'-matrix and det u = 2i. Another possible choice used here corresponds 
to uqr = uu = 1 and uir = uqi = defining C = TZ, the reactance matrix. The two 
representations are related as 

5 = (l + ^7^)(l-^7^)-^ (37) 

This identity holds for the exact matrices therefore it can be used as a check of the accuracy 
of the calculation by comparing the results using both schemes. 

At energies below the three-body breakup, the \l/c term is C^. Accordingly it can be 
represented by means of an expansion in the same C^ basis used for bound states, namely 

^, = Y,A„,,\mu) (38) 

mi/ 

From the above definitions we can construct the scattering state as 

* = ^ A^^\miy) + n^sT + ^^It (39) 

mv 

The solution of a scattering problem at a given energy requires the determination of 



the amphtude £ and t 
variational principle 



le linear coefficients Amu- To this aim we make use of the Kohn 



Iq that can be written as 



[C] = C--^{'^*\n-E\'^). (40) 

Clv^u til 
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The numerical implementation of the variational principle leads to a first order approxi- 
mation of the amplitude L obtained through the solution of a linear system of equations of 
size M + 1, where M is the size of the basis set for the expansion of the core part of the 
wavefunction. If we define an array of unknowns ({v4,„j,}, C) of dimension Af + 1, the linear 
system can be written as: 

Hm'u',mi/ Hjn'u',n^ \ I ^mu \ I —Hm'u,QP 

where H^' x stands for the matrix element 



(41) 



Hx',x = {x'*\n-E\x). (42) 

The second order estimate for C is then given by 

^2nd ^ ^ist ?_/^i«t*|^ _ ^1^1^*)^ (43) 

detw 

where \Ef^** is the wavefunction obtained solving the linear system of eq. pTl) . 

Let us now discuss in more detail the structure of eq. fHT]) . The top left part of the 
coefficient matrix, of dimension M x M contains the matrix elements used for the bound 
state calculation when the scattering state has the same quantum numbers as the bound state 
(compare it to eq. (I25p ). Otherwise specific states \mi>) having proper quantum numbers 
have to be constructed. The additional matrix elements needing to be computed are those 
between the C^ basis functions and the scattering functions, and among the scattering 
functions themselves, for a total of 2M + 4 different terms. The number of such extra terms 
grows linearly with the basis set size, and due to the functional form of f2^y and Qg^, they 
need to be calculated at every different choice of the system energy E. 

The application discussed above employ the HA basis in the expansion of the C^ \E'c term. 
Alternatively, \l'c could also have been expanded in terms of sole HH functions as 

^^ = J2^mki\mkl), (44) 

mkl 

where we have defined the ket 

\mkl) = L^^\(3p) exp [-/?p/2] ® \kl, OSTU). (45) 

After including a sufficient number of Laguerre polynomials, both expansions, in terms of HH 
or HA functions, are equivalent leading to the same value of C Example of this equivalence 
will be shown and discussed in the next Section. 

13 



B. Scattering below Break-up: Method 2 

An alternative approach considered is represented by a direct solution of eq. (fTO!) which 
represents a different form of the three-body Schroedinger equation. The bound state solu- 
tions have been discussed in Sect. II, and here we will discuss the scattering solutions below 
three-body breakup: E"^^ < E < 0. For this purpose it is important to determine the 
boundary conditions to be imposed to the functions {u^{p)}. Firstly, let us observe that at 
very large p the only open channel in the system of eqs. dlOl) is the lowest one, and that the 
system uncouples: 

-|-T, + Ui-E + Bn^ u, = 0. (46) 

At p = corresponds ui(0) = 0, whereas the boundary conditions at large p depend on the 
specific asymptotic forms of the hyperradial potentials Ui{p) and of the terms Bii{p). A 
detailed study of their asymptotic expressions will be object of a forthcoming publication 



ITJ I . For the purpose of this work it suffices to say that 

(-|^T, + U^-E + 5n) ^1 - (^ + ^2 + <P"'?j (p'/^«i), (47) 



where the wavenumber kp is defined from the relation: 

The boundary conditions associated with ui thus are 

/„x „ ,. / X ~ sin (A;„p) ^cos (kup) , ^. 

MiO =0, limnip^ui = \jf^+tan6 \jf^, 49 

all other u^ ^ sufficiently fast, as p ^ oo. Furthermore, the lowest adiabatic func- 
tion $5'"^^^(p, fi) -^ p'^/'^(l)d{r)\ST) at very large values of p [5|. Therefore, the asymptotic 
behavior of the scattering wave function in terms of the adiabatic basis results: 



^> = Y,nAp)^'r\p.n) ^ ur) 



!!!^iM + tan5^^^iM 



\ST). (50) 



P P 

In the limit p —>■ oo the relation ky-y ^ kpp holds as r is constrained by the finite size of the 
dimer wavefunction, therefore r/p <^ 1. Consequently eq. (!50l) represent the asymptotic 
limit of Qgrp + ta.n6Q^<^rp, for p —>■ oo. The full equivalence between the above expression 
for the asymptotic wavefunction and that one given by eqs. fl30|31l) can be established by 
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qR _ 


/ m 


V 2h^kp 


0^ - 

^''p,ST — 1 


1 m 


\j 2n?kp 



noticing that the ui constitutes the leading term in the expansion of Vt^rp and Vt^gj^ in terms 
of the small parameter r/p 6| , which yields 

(nf^|$i)-^^ + 0(p-^/^), (51) 

and 

(^frl*.) ^ ^^#^ + 0{p'){u > 1). (52) 

and a similar expansion yields for fi^j.. From the above discussion, we can define an alter- 
native asymptotic term $„ as combination of the following functions: 

;i-exp[-7p])''^^^$i(fi,p), (53) 

and 

l-exp[-7p])^^^*^$i(fi,p), (54) 

where the factor (1 — exp[— 7p])'' is introduced as usual to regularize the behavior of the 
functions for p ^ (in practical calculations we have set rj = 4), and the functions are 
normalized as in eq. (l32il . The same approach as in the previous Section can now be applied 
where. Accordingly the scattering wave function can be written as 

"^ = ^Bm^\miy) + n^ST + <S^lT (55) 

Mu 

where the asymptotic part is now given in terms of ^fsT ^^^ ^i st^ ^^^ ^^^ "^^^^ P^^^ ^c is 
expanded onto the HA basis times a set of C^ functions Uy{p). We can refer to this expansion 
as HA2. 

This approach is justified as the neglected terms in the r/p expansion of Vt^ and Vt^ do not 
carry flux and can be incorporated into the unknown term '^c- The approximated expression 
for the term ^a allows to speed up the calculation significantly as there is no need to 
calculate the overlap integrals between the HA basis functions and the asymptotic functions 
as in eq. fj^Tl) . On the other hand, its implementation suffers from the following problems. 
At intermediate distance the expansion on r/p of the asymptotic functions converges very 
slowly, resulting in a large number of HA functions which need to be taken into account. 
At large p, the implementation of the functions of eqs. fl30|3ip results in a very awkward 
behavior of the C^ term \E'c- Continuing the expansion of eq. (15T!) . for instance, it is possible 
to show that the next term is cos [/Cpp]/p'^/^, which imposes the asymptotic behavior that 
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the function ui has to reproduce. This particular functional form is very slow decaying, and 
it is particularly hard to reproduce with a polynomial expansion. This problem is further 
enhanced by the presence of oscillations associated with cosine and sine terms. 

In order to solve the linear system fHT]) taking into account the oscillatory behavior 
of the hyperradial functions for large p values, we have implemented a Discrete Variable 



Representation (DVR) scheme [18|] rather than the standard variational approach. In a 
previous work [l9| we have shown how to combine the variational Kohn principle with a DVR 
scheme, for the case of a two-body system, which corresponds to a single one-dimensional 
differential equation. In this work we have a set of Na one-dimensional coupled differential 
equations. Therefore we define a {N^M + 1) x [N^M + 1) unitary transformation matrix 
U which is a direct product of N^ + 1 matrices 

U = U^'^®U^'^®...®1, (56) 

where lA^'^ is a M x M unitary matrix associated to a customary one- dimensional DVR of 
size NovR = M built in p: 

W^ = Lf(t,)exp[-t,/2]V^, (57) 



where tj and Wj are the appropriate quadrature points and weights. By mean of a parameter 
/?, the end quadrature point t^^^^ can be associated to different physical values Pmax, by 
setting tj = (3pj. In this fashion we can constrain the quadrature points to be distributed 
between and pmax- 

IV. NUMERICAL APPLICATIONS 

In order to illustrate the method outlined in the previous Sections we present two ap- 
plications to the n — d system in a quartet state (5* = 3/2). The potential energy of the 
system is taken as the sum of three pairwise potentials. We consider the MT-III interaction 
Vmt-iii for which benchmarks results exist in the literature [9|. It reads: 

VMT-iii{r) = (1438.72 exp [-3.11 r] -626.885 exp [-1.55 r])/r. (58) 



E 



To make contact with the results of Ref . [6| , we have also used the Gaussian potential (named 

Vg): 

Vcir) = -66.327 exp [-(0.64041 r)^], (59) 
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For both potentials we assume nuclear distances in fm and energies in MeV. The nucleon 
mass used is such that fv^ /m = 41.47 MeV fm^. Furthermore, we consider both potentials 
as acting only on the / = two-body partial wave. 

The potential Vq supports one deuteron bound state, with zero angular momentum, of 
energy E2b = —2.22448 MeV. The zero-energy scattering length is a^ = 5.4208 fm, whereas 
for the MT-III potential the values are E2b = -2.23069 MeV, and a^ = 5.5132 fm. 

For the potential Vq we consider the three-body system with quantum numbers 11 = a, 
T = 1/2 and S = 1/2, whereas for Vmt-iii H = a, T = 1/2 and S = 3/2. As the potentials 
are projectors on s— wave, the index / in eq. ( TT^ is restricted to the value / = 0, and the 
index k can take the values /c = 0, 2, 3, 4, 5, . . . , oo in the first case and k = 1,2,3,4,5, ... ,oo 
in the second case. 

A. Bound states 




FIG. 1: The top panel shows the lowest adiabatic curves Uy{p) for different values of p. In order 
to display the behavior at large p the curves are multiplied by a factor p^. The lowest curve 
thus tends to E2bP^, and the others to the spectrum h?K{K + 4)/m, with K = 0,2,3, 4, .. . for 
u = 2,3,4,5, .. .. The bottom panel shows the convergence of the lowest adiabatic curve Ui{p) as a 
function of the number of HH used in the expansion of eq. (|13p . The asymptote at i?2b = —2.2245 
MeV is plotted for comparison. 



Figure [T] shows, in the upper panel, the lowest adiabatic curves Uy{p) calculated for 
the Vg potential. In order to highlight their asymptotic behavior, the curves have been 
multiplied by a factor p^ . The lowest curve f/i(p) thus tends to the deuteron energy times 
p^ , whereas the upper curves tend to the free HH spectrum, that is 4k{k + 2)h'^/m with 
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k = 0, 2, 3, 4, 5, . . . for z/ = 2, 3, 4, 5, 6, . . .. The value A; = 1 is not allowed as there is no 
completely symmetric HH with k = 1 and I = 0. Subsequently, the adiabatic function 
$i(p, i7) tends to the deuteron wavefunction, whereas ^^, z/ > 1, to the HH functions, with 
the appropriate normalization factors. The lower panel shows the convergence of the lowest 
curve Ui{p) as a function of the number of HHs employed in the expansion of eq. ( 1131) . It 
shows that the larger p becomes, the larger the expansion basis must be in order to properly 
describe the function $i. In practice, the radius of convergence of expansion flT3l) increases 
rather slowly when the basis set size is increased. The reason for this behavior is that 
when p is increased the function $i becomes more and more localized in the hyperangular 
phase-space, therefore its description by means of the HH requires a larger and larger basis 
set size. This behavior is not connected with any particular feature of the potential used 
in this specific calculation but it can be considered a general one, as it is induced by the 
geometric localization of the deuteron wavefunction in connection with the HH expansion. 
The thick curve is the solution of eq. ( l2TT) starting at p = 20 fm. For large values of p, the 
corresponding eigenvalue reproduces the two-body binding energy E2b- 

The description of a three-nucleon bound state using a central potential has to been 
taken as a homework problem and preliminary to check the usefulness of the HA basis to 
treat scattering states, in comparison to the HH expansion. The Vg potential predicts two 
bound states in the three-body system, a very deep ground state and a very shallow excited 
state. Table [I reports the convergence patterns for the upper bounds E^ and E2 to the two 
bound states supported by the potential Vg, as a function of the number A^ of HA and HH 
basis elements. The HA functions were expanded in 80 HH functions which is the number 
required for the HH expansion to describe accurately the deep and shallow bound states. 
The number of HH functions necessary to obtain a full convergence of the energy for the 
deep bound state is much smaller, around 10 functions. The most striking feature to be 
observed in the table is the much rapid convergence of the HA basis expansion compared 
to the HH. Not only full convergence can be achieved with a basis which is one order of 
magnitude smaller, but already the inclusion of only one basis element yields an energy for 
the excited state within 90% of its converged value. 
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n=l n=2 

A^ HH HA N HH HA 

1 -21.5808 -22.0520 1 0.0620 -2.3484 

2 -21.9567 -22.0850 4 -0.9576 -2.3627 

3 -22.0694 -22.0873 10 -2.0348 -2.3632 

4 -22.0805 -22.0874 20 -2.3036 -2.3632 

5 -22.0852 -22.0874 30 -2.3474 -2.3632 

6 -22.0869 -22.0874 40 -2.3582 -2.3632 

7 -22.0872 -22.0874 50 -2.3615 -2.3632 

8 -22.0873 -22.0874 60 -2.3626 -2.3632 

9 -22.0874 -22.0874 70 -2.3631 -2.3632 

10 -22.0874 -22.0874 80 -2.3632 -2.3632 

TABLE I: Patterns of convergence for the three-nucleon bound states obtained with the Vq po- 
tential, as a function of the number A^ of hyperangular basis functions included in the expansion. 
The HA basis elements were calculated with 80 HH, (3 = 1.6 fm~^, and 33 Laguerre polynomials 
were employed in the expansion of eq. (j24p . Note the different scales for the ground and excited 
state patterns of convergence. 

B. Scattering States 

In the following, results obtained combining the HA basis expansion with the expressions 
of eqs.( !30|3Tl) are given and will be referred to as HAl. Table HTl reports the full patterns 
of convergence of the L = 0, S = 3/2 MT-III phase shift 6, at Ecm = 1 MeV, as a function 
of the number of Laguerre polynomials Np used in expanding the hyperradial functions in 
eq. fl38l) and the number Na of adiabatic channels included. The HA functions have been 
calculated using 200 HH functions. This number of HH functions is sufficient to accurately 
describe the phase shifts below the three-body breakup. From the table it can be seen that 
the convergence requires a rather high number of HA basis elements, more than 100, whereas 
12 Laguerre polynomials are enough to achieve final convergence. 

In order to analyze deeply the pattern of convergence, in Table IIIII results obtained by 



means of the HH expansion 20] are compared to those obtained with the HA approach. In 
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Np\Na 20 40 60 80 120 160 200 

5 -55.974 -55.912 -55.902 -55.898 -55.897 -55.896 -55.896 

9 -55.937 -55.879 -55.870 -55.867 -55.865 -55.864 -55.864 

13 -55.932 -55.878 -55.868 -55.865 -55.864 -55.863 -55.863 

17 -55.934 -55.878 -55.868 -55.865 -55.863 -55.863 -55.863 

21 -55.932 -55.878 -55.868 -55.865 -55.864 -55.863 -55.863 

25 -55.933 -55.878 -55.868 -55.865 -55.864 -55.863 -55.863 

29 -55.932 -55.878 -55.868 -55.865 -55.864 -55.863 -55.863 

33 -55.931 -55.878 -55.868 -55.865 -55.864 -55.863 -55.863 

TABLE II: Convergence of the phase-shift 6 in function of the number of Laguerre polynomials 
Np (see eq. (j24p ) and of the size Na of the HA basis set, at an incident energy of -E = 1.00 MeV. 
The HA basis is calculated with 200 HH elements. The non-linear parameter was fixed to /3 = 1.9 
fm~^. 

each row of the table N^ indicates the number of HH functions used in the calculation and 
the number of HA functions used calculated using 200 HH functions. As already pointed out, 
for the special case of Nh = Na the two expansions are equivalent and the results become 
identical, provided that a sufficiently high number of Laguerre polynomials is employed to 
describe the {u^ip)} set of functions. Therefore the equivalence can be seen in the last row of 
the table in correspondence with Na = 200 (in some cases the equivalence is reached already 
at A^^ = 160). For the case oi E = 2.00 MeV, two patterns of convergence are shown for two 
different HA bases, obtained with 120 HH and 200 HH, respectively. Here the equivalence 
can be seen also at Na = 120. In this energy range there is little difference in the results 
obtained with the two bases, for example when 20 or 40 HA basis elements are employed. To 
be noticed that the results shown in Table UTTl present a different pattern of convergence with 



respect to the ones given in Table 2 of Ref 



20| : the reason is that in the previous paper the 



5— matrix representation was chosen for the matrix u, whereas in this work the i?— matrix 
was preferred. The two choices are equivalent and lead, once convergence is achieved, to the 
same results. We can conclude that although there is some improvement, the table shows 
that the convergence is not speed up significantly by transforming the HH basis into the 
HA basis. This suggests that the HA basis does not provide as an optimized basis for the 
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scattering problem as it does for the bound state problem. 





0.20 MeV 


1.00 MeV 


2.00 MeV 


Na 


HH 


HAl 


HH 


HAl 


HH 


HAl 




120 200 


20 


-28.263 


-28.312 


-56.913 


-55.931 


-70.741 


-71.594 -71.597 


40 


-28.201 


-28.299 


-55.948 


-55.878 


-71.701 


-71.501 -71.500 


60 


-28.306 


-28.295 


-55.922 


-55.868 


-71.508 


-71.485 -71.483 


80 


-28.296 


-28.294 


-55.872 


-55.865 


-71.483 


-71.480 -71.478 


120 


-28.294 


-28.294 


-55.865 


-55.864 


-71.476 


-71.476 -71.475 


160 


-28.294 


-28.294 


-55.863 


-55.863 


-71.474 


-71.474 


200 


-28.294 


-28.294 


-55.863 


-55.863 


-71.474 


-71.474 



TABLE III: Convergence of the phase-shift 6 at three different energies below break-up threshold 
for the MT-III potential, in function of the size A^ of the basis. The patterns of convergence for 
the HH and HAl methods are shown for comparison. The HA basis was calculated employing 200 
HH basis elements. For E = 2.00 MeV the calculation with 120 HH basis elements is also shown. 
All calculations employed 33 Laguerre polynomials (see Table [H]), and (3 = 1.9 fm~^. 

Table [IV] shows the convergence pattern for the phase-shift at E=1.00 MeV, obtained 
using the HA2 expansion for the asymptotic term. As anticipated in the previous Section, 
in order to obtain stability in the phase shift, we have employed a much larger and finer 
hyperradial grid, consisting of 4153 points, distributed up to p = 2000 fm. At the same 
time the HA basis set and associated eigenvalues were obtained with a bigger number, up 
to 2000, of HH basis functions, or by solving the asymptotic differential eq.( l2T1) for p > po 
(po = 40 fm). This calculation has been performed using the Laguerre polynomials as an 
expansion basis for the hyperradial functions. As anticipated, the polynomials are not 
an appropriate choice to reproduce the long range oscillatory behavior of the hyperradial 
functions. This can be seen from the poor convergence pattern in terms of Np as the 
number of HA functions increases. For A^^ > 8 more than 100 polynomials are necessary. 
Furthermore, the convergence pattern is also poor relative to the increase of the number of 
HA basis elements. Differences with results of Table HT] are remarkable. 
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Np\NA 
21 
41 
61 
81 
101 
121 



16 



24 



32 



36 



40 



-57.753 -57.231 -57.063 -57.037 -57.230 -57.028 -57.027 
-57.638 -56.915 -56.581 -56.511 -56.489 -56.484 -56.480 
-57.628 -56.868 -56.456 -56.348 -56.310 -56.300 -56.293 
-57.627 -56.858 -56.414 -56.281 -56.228 -56.214 -56.204 
-57.626 -56.855 -56.399 -56.251 -56.188 -56.169 -56.156 
-57.626 -56.853 -56.393 -56.237 -56.166 -56.145 -56.129 



TABLE IV: Convergence of the phase-shift 5, using the HA2 method, in function of the number 
of Laguerre polynomials Np (see eq. (I24p ) and of the size Na of the HA basis set, at an incident 
energy of E = 1.00 MeV. The HA basis is calculated with 2000 HH elements. The non-linear 
parameter was fixed to /? = 1.9 fm^^. 



FIG. 2: The phase-shift 6 in terms of different choices of the non-linear parameter f3 and of the 
size of the expansion in Laguerre polynomials. The top panel shows the convergence for expansion 
HAl, and the bottom panel for expansion HA2. Note the different scales on the y— axis of the two 
graphs. 

Figure [2] shows the effect on the phase-shift of varying the non-linear parameter jS. The 
upper panel shows results for the HAl expansion, whereas the lower panel refers to the HA2 
expansion. Different sizes of the Laguerre basis are shown. In principle, for a complete basis 
set, that is A'"p = oo, there should be no effect in varying the parameter /3. When the basis 
set is finite, the stability of the result, in this case the phase-shift, with respects to changes 
of /3 is a measure of the completeness of the expansion. In particular, by comparing the 
upper and lower panels, one can see that the HAl polynomial expansion of the functions 
Uu{p) is much more effective than for the case HA2 (also note the different scales of the 
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Pma.y\NDVR 


100 


150 


200 


250 


300 


350 




-ist 


cynd 


200 


-56.179 


-56.161 


-56.159 


-56.159 


-56.159 


-56.161 


-56.160 


400 


-56.124 


-56.100 


-56.095 


-56.093 


-56.092 


-56.091 


-56.092 


600 


-56.096 


-56.089 


-56.085 


-56.084 


-56.084 


-56.085 


-56.083 


800 


-56.119 


-56.089 


-56.084 


-56.083 


-56.082 


-56.080 


-56.081 


1000 


-56.162 


-56.087 


-56.082 


-56.081 


-56.081 


-56.082 


-56.081 


1200 


-56.149 


-56.088 


-56.082 


-56.081 


-56.081 


-56.082 


-56.081 


1400 


-56.106 


-56.084 


-56.082 


-56.081 


-56.081 


-56.077 


-56.080 


1600 


-56.154 


-56.090 


-56.082 


-56.081 


-56.081 


-56.082 


-56.080 



TABLE V: Convergence of the phase-shift 5 sX E = 1.00 MeV, using the HA2 method, for the 
MT-III potential, as a function of the number N^VR of DVR points employed, and of the last grid 
point Pmax- Convergence is shown for the second order estimate of S for all values of M, but the 
last, where both first and second order are shown. 

J/— axis). In the first case, the expansion with 17 polynomials is completely unaffected by 
changes in /?, whereas in the second case even a basis set as large as 120 polynomials yields 
significantly different results with different choices of /3, indicating that the result is far from 
convergence. 

In order to circumvent this problem we use the DVR technique in the hyperradius variable. 
Table IVl shows the convergence, in terms of different choices of pmax and the number of DVR 
points employed, of a case calculation, with 40 adiabatic functions, for the MT-III potential 
and E = 1.00 MeV. For the biggest case {Ndvr = 350), we show both the first and second 
order values of the phase-shift obtained by using the Kohn Variational Principle. In order 
to obtain a good convergence of the second order value it is important that the integral in 
eq. fj43l) is calculated with a very high numerical accuracy. The hyperradial grid used in the 
calculation consists in more than 4000 grid points up to p = 2000 fm. The use of the DVR 
technique allowed for stable results in terms of the hyperradial expansion. The use of 350 
DVR points is equivalent to a calculation with 350 Laguerre polynomials which in general 
is much more involved to be carried. However the number A^^ = 40 of HA functions used 
in this calculation is not enough to well describe the phase shift. At E=1.00 MeV the HAl 
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method as well as the HH method predict S = —55.863 degrees to be compared to the result 
of the HA2 method, 6 = —56.081 degrees, using Na = 40. In order to have an stable result 
for 6 using the HA2 method, the value Na = 120 has to be considered and Ndvr > 350 since 
the number of DVR points has to be increased as A^^ increases. The dimension of the HA2 
problem is A^^ x N^vr and is clear that very soon the problem becomes computationally 
unsustainable, unless exceptional computational resources are considered. 





0.20 MeV 


1.00 MeV 


2.00 MeV 


Na 


HAl 


HA2 


HAl 


HA2 


HAl 


HA2 


4 


-28.364 


-29.065 


-56.136 


-57.625 


-72.344 


-71.988 


8 


-28.340 


-28.739 


-56.038 


-56.852 


-71.965 


-71.871 


12 


-28.328 


-28.604 


-55.984 


-56.545 


-71.770 


-71.437 


16 


-28.319 


-28.532 


-55.947 


-56.385 


-71.660 


-71.210 


20 


-28.312 


-28.487 


-55.922 


-56.286 


-71.597 


-71.070 


24 


-28.308 


-28.456 


-55.906 


-56.218 


-71.558 


-71.975 


28 


-28.304 


-28.434 


-55.895 


-56.169 


-71.534 


-71.907 


32 


-28.302 


-28.417 


-55.888 


-56.133 


-71.518 


-71.855 


36 


-28.300 


-28.404 


-55.882 


-56.104 


-71.507 


-71.815 


40 


-28.299 


-28.394 


-55.877 


-56.081 


-71.500 


-71.783 


Table [TTT 


-28.294 


-55.863 


-71.474 



TABLE VI: Patterns of convergence for the two different choice of the asymptotic term, in terms 
of the number Na of HA basis elements, at three different energies. The MT-HI potential has been 
used. The last row reports the converged values from Table Hill The columns refers to a choice of 
P = 1.9 fm^^ for HAl, and pmax = 1200 fm for the HA2 expansion. Moreover, the HAl values 
are associated to a calculation with 200 HH, whereas the HA2 to a calculation with 2000 HH. The 
HA2 results have been obtained with the DVR scheme. 



Table |Vl] compares the convergence patterns for 6 at three different energies for the two 
suggested choices for the asymptotic term \E'a, namely the one in eqs. fl30|31l) . referred to 
as HAl, and the one in eqs. fl53ll54p . referred to as HA2, in terms of the number A^^ of 
adiabatic channels. Due to the very large basis sets required to obtain convergence with 
the HA2 term, the pattern of convergence is limited to few channels, less than required to 
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obtain a full convergence. The last row reports the converged values from Table lllli It is 
possible to see that the expansion HAl converges faster towards the final number, whereas 
expansion HA2 moves rather slowly. The reason is the difference in the treatment of the 
asymptotic wavefunction. In the HAl method, as well as in the HH method, the asymptotic 
configuration described by \E'a is reached at intermediate distances. Conversely, in the IIA2 
the the configuration described by \l/a is reached at much larger values of p. Furthermore, 
at intermediate distances, in order to reproduce the correct behavior a big number of HA 
functions are needed. 

The following figures present important characteristics of the hyperradial functions used 
in the expansion HA2. 





FIG. 3: The functions Uy{p) {u = 1,2,3,4) a± E = 2.00 MeV. The top panel shows the short- 
range region. Some of the functions were magnified by the factor shown in the legend. The circles 
represent the DVR amphtudes at the DVR grid points. The lower panel shows the long range 
region. In order to highhght the asymptotic behavior of each function, ui is multipUed by p^'^, 
and U2, U3 and U4 by p^. In the bottom panel ui is also magnified by a factor 50000. 



Fig. [3] shows the functions u^{p) calculated with A^^ = 4, N^vr = 300 and pmax = 1200 
fm. The dots indicate the DVR amplitudes at the DVR points, whereas the lines represent 
the u^{p) functions obtained by back-transforming to the original polynomial basis. The top 
panel displays the short-range region (0 < p < 20 fm), where the function ui is predominant. 
The bottom panel shows a part of the long-range region (100 < p < 300 fm). Here the 
situation is drastically different, and the functions U2, M3 and M4 have a much larger amplitude 
than Ui (which is magnified by a factor 50000). Also, in order to highlight the asymptotic 
behavior, ui is multiplied by p"^/^, and U2, u^ and M4 by p^. The most striking feature are 
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the oscillations present in all curves. This behavior is a consequence of the decomposition of 
the asymptotic configuration in terms of HA functions. This resulting peculiar long range 
behavior is the cause of the very slow convergence of the phase-shift shown in Tables IIV[ FV] 
IVll The behavior obtained for the curves u^ is the one expected by the analytical expansion 
of the asymptotic terms indicated in eqs. fl51ll52p . 




FIG. 4: The function Uii{p) (see eq.®) obtained with the method HAl (continuum line) and 
HA2 (dotted line). In the first case the function u^ is short range, decaying exponentially with p, 
whereas in the second case it shows the long range oscillations. 

Figure m compares the hyperradial function Ui{p) obtained with method HAl and HA2. 
In particular it highlights as the former is short range and exponentially decaying with p, 
compared to the latter which is oscillating as indicated in eq. (ISTl) . 



As mentioned in the Introduction, in Ref. 6|] the phase shift for the potential Vq has 
been calculated from eq. fITOl) in the so-called uncoupled adiabatic approximation (UUA) 
retaining one hyperradial function. Namely, the following equation has been solved: 



2m 



ui{p)=0 



(60) 



with the asymptotic condition Ui{p) -^ sin{kp + 6 + 37r/2) as p ^ oo. Besides the factor 
37r/2, this is equivalent to the method HA2 given in the previous section taking into account 
one HA function. 

In Figure [5] we show the phase-shift S{E). The dots represent fully converged results 
obtained with the HAl expansion, whereas the continuum line represent results obtained by 
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FIG. 5: Elastic deuteron-nucleon phase-shift below three-body break-up (marked by the dotted 
line) for the Vq potential. The full line corresponds to a calculation retaining one HA basis element, 
and using the HA2 method. The dots correspond to the full calculation. 

including just one adiabatic function in the expansion HA2. It is possible to notice that the 
UUA provides a very good first order estimate of the phase-shift. However, the deviation 
from the complete expansion can be as big as 10%. Also notice that in Figure [5] the phase- 



shifts have been normalized so that 6{E = 0) — 6{E = oo) 
trimer states. 



360, as there are two bound 



V. CONCLUSIONS 



In this paper we have investigated the capability of the HA basis to describe scattering 
states in a three-nucleon problem. The basis was generated from the hyperangular Hamil- 
tonian by means of an expansion in HH functions. We have shown the complete equivalence 
between the adiabatic basis generated using A^ HH functions and the HH basis of dimension 
A^. This equivalence provides a useful benchmark when the convergence of the quantities of 
interest is studied in terms of the number A^^ of adiabatic functions. For example, for bound 
states it is well known that A^i << N suffices for the convergence of the binding energies. 
One goal of this paper was to investigate whether the same relation holds for scattering 
states. In particular, we studied the convergence of the L = phase shift S corresponding 
to a process in which a nucleon collides a deuteron at low energies in the state S = 3/2. For 
this purpose we have used the MT-III potential. 

In the calculation of the phase shift using the HA basis we have followed two different 
procedures. They were both based on a decomposition of the scattering wavefunction as 
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a sum of two terms. One term describes the configurations when the three particles are 
all close to each other and goes to zero as the interparticle distances increase. The second 
term describes the asymptotic configurations and has been regularized so that goes to zero 
as y —>■ 0. In the first procedure the HA basis has been used to expand the short range 
part of the scattering wave function. The second order estimate of the phase-shift has been 
obtained from the Kohn variational principle. A similar approach has been used before 
with the HH basis. Therefore, a detailed comparative analysis of the convergence patterns 
was possible. The conclusion is that the number of basis elements needed to achieve a 
comparable level of convergence for the phase-shift is of the same order for the two bases, 
that is Na ^ N, which is a surprising difference with respect to what happens in bound state 
calculations. A possible explanation could be the following. In bound state calculations the 
wavefunction expansion benefits from the initial optimization performed by constructing the 
HA basis. Conversely, in scattering state calculations the solution of the linear system of 
eq.( l4T]) requires a different short range behavior in the HA basis elements due to the presence 
of the terms fl^gj^ and Qgj^ in the short distance region. 

The second procedure considered was based in a direct solution of the system of equations 
for the hyperradial functions given in eq. flTOj) . This method however suffers from the follow- 
ing complications. The hyperradial boundary conditions to be imposed are those required 
to reconstruct the asymptotic configuration given by the functions defined in eqs. (13011311) . 
For very large values of p the boundary conditions are simple and are given by eq. (l49l) 
for the lowest function {ly = 1). All other functions go to zero as p ^ oo. This means 
that the solution of the linear system has to be obtained over a very extended hyperradial 
grid. Moreover, the adiabatic potentials and functions have to be accurately known in the 
grid. In the present work we have solved partially the numerical difficulties associated to 
the solution of eq. (TTOl) introducing a variational DVR procedure. 

From the present study we can conclude that the use of the HA basis in the description 
of scattering states is not as advantageous as for bound states. The main drawback is that 
then number of basis elements required to reach convergence is not as low (in proportion) 
as in bound state calculations. Secondly, a number of numerical problems arise from the 
need of calculating the adiabatic curves and the associated basis elements at large distances. 
Further studies to improve the description of scattering states using the HA expansion are 
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at present underway. 
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